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After deriving the isothermal Hellmann-Feynman theorem (IHFT) that is suitable for mixed states in a canon- 
ical ensemble, we use this theorem to obtain the isothermal magnetic-field sweep theorems for the free, aver- 
age and trapping energies, and for the entropy, specific heat, pressure and atomic compressibility of strongly- 
correlated ultra-cold quantum gases. In particular, we apply the sweep theorems to two-component Fermi gases 
in the weakly-interacting BCS and BEC limits, showing that the temperature dependence of the contact parame- 
ter can be determined by the variation of either the entropy or specific heat with respect to the scattering length. 
We also use the IHFT to obtain the Virial theorem in a canonical ensemble, and discuss its implications for 
quantum gases. 
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I. INTRODUCTION 

Dilute samples of ultra-cold quantum gases are well- 
described by zero range interactions, a consequence of which 
is the existence of some universal relations that govern the 
behavior of these systems [1]. The first examples of such 
relations were originally derived by solving the many-body 
Schrodinger equation flj-Hl, and they relate energy and pres- 
sure of the system to a single parameter that comes from the 
short-range physics. What makes these relations remarkable 
is that they hold for any finite-energy state of the system, it 
does not matter whether the system is few- or many-body, su- 
perfluid or normal, weakly or strongly interacting, in equilib- 
rium or nonequilibrium, at zero or finite temperature. These 
universal relations have more recently been rederived using 
many other approaches including the quantum field theoretical 
techniques lll|5|-|9t], and have also been verified via numerical 
few-body calculations ifloll . 

There is one common element in these universal relations: 
they all involve a single parameter called the contact. In the 
case of two-component Fermi gases, the contact parameter is 
given by the coefficient of the \/k A tail of the momentum dis- 
tribution, and therefore, it measures the number of fermion 
pairs with small separations. It has been found that the con- 
tact parameter also appears in many other short-range (high- 
momentum) or short-time (high-frequency) properties of the 
system ill lWlfJl . This parameter has recently been measured 
in an ultracold 40 K gas via the measurements of the high- 
momentum tail of the momentum distribution and of the high- 
frequency tail of the radio-frequency signal fl7ll . and also 
measured in an ultracold 6 Li gas via the measurements of the 
static structure factor Hill . The measured temperature and 
scattering length dependence of the contact parameter com- 
pare well with the theoretical predictions 1 19112011 . 

In particular, the so-called adibatic sweep theorem for the 
energy [3], which has recently been experimentally veri- 
fied fl7tl . relates the contact to the change in total energy of the 
system when the atom-atom scattering length (a) is changed 
adiabatically, i.e. with zero heat transfer. In atomic systems, 
the value of a can be tuned at will from very large and nega- 
tive to very large and positive values, thanks to the presence 
of magnetically induced Feshbach resonances, and an adia- 



batic (constant entropy) sweep is accomplished by changing 
the bias magnetic field over a time scale longer than the relax- 
ation time 11711 . 

Motivated by these earlier works, here we generalize a 
number of these universal relations to finite temperatures, and 
the rest of the manuscript is organized as follows. In Sec. [II] 
we derive the isothermal Hellmann-Feynman theorem (IHFT) 
that is suitable for mixed states in a canonical ensemble. Then, 
we use this theorem in Sec.[III]to obtain the isothermal (con- 
stant temperature) sweep theorems (1ST) for the free, average 
and trapping energies, and for the entropy, specific heat, pres- 
sure and atomic compressibility of strongly-correlated ultra- 
cold Fermi gases. We also use the IHFT to derive the Virial 
theorem in a canonical ensemble in Sec.[lV] A brief summary 
of our conclusions is given in Sec. [VI 



II. HELLMANN-FEYNMAN THEOREMS 

For our purpose, we need to derive the Hellmann-Feynman 
theorem (HFT) suitable for a canonical ensemble of a fixed 
number N of identical particles that are in thermal equilibrium 
with a heat reservoir, which is discussed next. 



A. HFT for pure states 

For any stationary normalized eigenvector \ip n (X)) with the 
corresponding eigenvalue E n (X) of a Hamiltonian H(X), i.e. 
(VV1IVV1) = 1 and E n = (ip n \H\ip n ), the usual HFT states 
that 



dEn 
dX 



dH 



ax 



(1) 



where A is an arbitrary real parameter, explicitly appearing in 
H. This well-known theorem has proved to be very useful in 
many fields ranging from quantum chemistry, quantum statis- 
tics and many-body physics to molecular physics, with a lot 
of applications. 

In particular, the HFT has recently been used to derive a 
number of exact relations for strongly-correlated systems with 



2 



short-range interations, in the context of ultra-cold quantum 
gases ||5l |g, 12611 . Since this theorem applies only for pure 
states, these exact relations strictly hold at zero temperature 
for any change or equivalently at finite temperatures for an 
adiabatic change. In order to generalize these relations to fi- 
nite temperatures, here we derive the IHFT for a canonical 
ensemble of a fixed number of identical particles that are in 
thermal equilibrium with a heat reservoir. 

B. IHFT for mixed states in a canonical ensemble 



for strongly-correlated ultra-cold quantum gases, which is the 
main purpose of this manuscript. 



HI. ISOTHERMAL SWEEP THEOREMS 

In this section, we use the IHFT given in Eq. (f3j) to derive 
the isothermal sweep theorems (1ST) for the free, average and 
trapping energies, and for the entropy, specific heat, pressure 
and atomic compressibility in a canonical ensemble. 



For this purpose, we assume the temperature T is fixed, 
and the Hamiltonian H is a function of some arbitrary real 
parameter A. There are two complimentary ways to obtain the 
isothermal HFT in a canonical ensemble. 

The simplest way follows from the definition of the statis- 
tical Helmholtz Free energy F, i.e. F = — T In Z in a unit of 
ks = 1 here and throughout, where Z = Tr(e~ H / T ) is the 
canonical partition function of the system. Then, it is easy to 
show that (see also JHHHl) 



dF _ /dH\ 
SA ~ \BX/ 



(2) 



where the ensemble average of any operator is defined by 
(A) = Tr(pA) and p = e~ H / T /Z is the density operator 
in the canonical ensemble. We recall that the derivatives with 
respect to A are evaluated at fixed T here and throughout. In 
addition, using the definitions of the thermodynamic Free en- 
ergy, average energy and entropy, i.e. F = E — TS, E = (H) 
and S = — Tr(plnp) = —dF/dT, respectively, Eq. (0 can 
be written as 
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This is the isothermal HFT for mixed states in a canonical 
ensemble, and it will play a central role in the remaining parts 
of this manuscript. 

An alternative way of deriving Eq. OJ is as follows. We 
start directly from the definition of E, i.e. E = Tr(pH), and 
evaluate the derivative to obtain 
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In the intermediate steps, we used Tr(<9e H / T /dX) — 
-(l/T)Tr(e- H / T dH/d\) and Tr(Hde~ H / T /dX) = 
-(l/T)Tr(He- H / T dH/dX), both of which can be derived 
by writing Tr(^4) = J^ni^l^l^n) where \ip n ) is the normal- 
ized eigenvector and E n is the corresponding eigenvalue of H. 
In addition, we used de- E -l T jdX = {ip n \de- H / T / dX\ij) n ) , 
which can be obtained from the usual HFT given in Eq. (Q]). 
For fixed T, starting from the definition (A) = Tr(pA), it is 
easy to show for the observables that commute with H, i.e. 
[H, A] = 0, that d(A)/8T = ({H - E)A)/T 2 . Therefore, 
Eqs. (0 and © are equivalent. 

Having derived the IHFT for a canonical ensemble, next 
we use this theorem to obtain a number of exact relations 



A. 1ST for the free, average and trapping energies 

As a first application of the IHFT, we consider the Hamilto- 
nian that describes two-component Fermi gases in an external 
potential U a (r), i.e. 



ft 2 V 2 
2m„ 



+ U a (r) 



Vv(r) 



g / d 3 r^(r)^(r)^(r)Vt(r), 



(5) 



where ip^(r) creates a pseudospin-er fermion with mass m a , 
and g > is the strength of the short-range interaction. In 
the following discussions, the numbers of f and J, fermions 
need not be equal. As usual the theoretical parameter g can be 
written in terms of the experimentally more relevant scatter- 
ing length a via l/g = -MV/{4Trh 2 a) + (M/h 2 )Y /k (^/k 2 ), 
where M — 2ra^m^j (m^ + m^) is twice the reduced mass 
of j" and | fermions, and V is the volume. This equation 
gives g = —4:n 2 h 2 a/(MVT: — 2MVak c ), where fc c is the 
momentum-space cut-off used in order to evaluate the ultravi- 
olet divergent sum over k. 

Following the recent work Jst], the contact parameter C can 
be defined as 



dH \ 
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3 r(^|(r)^(r)^(r)^t(r)), (6) 



(7) 



where we used dg/da = —MVg 2 /(4irh 2 a 2 ). Note that C is 
an extensive quantity that depends on both a and T. This def- 
inition guarantees B5D that the average energy E of the system 
is of the desired form J2,[Ml 5 



E 



h 2 k 2 
2m n 



n CT (k) 



h 2 c 

AnMa ' 



(8) 



where n a (k) is the momentum distribution of a fermions. As 
emphasized in [2], this relation holds for any finite-energy 
state of the system, it does not matter whether the system is 
few- or many-body, superfluid or normal, weakly or strongly 
interacting, in equilibrium or nonequilibrium, at zero or fi- 
nite temperature. Combining Eq. (O with the IHFT given in 
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Eqs. (O and (01, we obtain 

dF h 2 C 
da 
dE 
da 



AirMa 2 

H 2 
A-nMa 2 



6 1 dT 



(9) 
(10) 



Here, the derivatives with respect to a are evaluated at fixed 
T, and therefore, Eqs. (0 and ( fTOb correspond to the 1ST for 
the free and average energies, respectively. Compared to the 
zero temperature expression for any change or equivalently 
the finite temperature expression for an adiabatic change, 
i.e. the adiabatic sweep theorem for the energy |Q1 [H, the 
main difference in Eq. ( fTOb is an extra TdC/dT term. Sim- 
ilarly, we can also calculate dF/dk c — h 2 C / (2irM) and 
dE/dkc = [h 2 /(2nM)}(C - TdC/dT), where we used 
dg/dk c = -MVg 2 /(2Tr 2 h 2 ). 

For homogenous systems (no external potential), the T de- 
pendence of C has been recently calculated in the low, inter- 
mediate and high T regimes as [19] 



C = C Q 
C = Co 

T ' 



C 



ai T\ (T«T C <7» 
a 2 T 2 , (T c < T < Tf) 

(T»{T F ,T a }) 



where, up to the leading orders in a, Cq 

9V^7r 6 n 2 a 2 /(40T|) and a 2 = 



oei 

l)Tr 3 k F ri 2 a?/(5T' F ) in the BCS limit; C* 
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(ID 

(12) 
(13) 

= 47r 2 n 2 a 2 , 
-4(7 In 2 - 
2imja and 



2ir 6 ^Jirkp / a(da m /da)n/T F in the BEC limit; and 
8w 2 n 2 /M for all couplings. Here, kp is the Fermi mo- 
mentum, T c is the critical temperature for superfluidity, Tp — 
h 2 k 2 F /{2M) is the Fermi temperature, T a = h 2 /(Ma 2 ) is 
the binding energy of two fermions, n = N/V = fc F /(37r 2 ) 
is the density of fermions, a m is the dimer-dimer scattering 
length between Cooper molecules (e.g. a m — 0.6a when 
m<|- = mi), and the weakly-interacting BCS and BEC lim- 
its are characterized by kpa — > 0~ and kpa — > + , respec- 
tively. Note that our definition of C is larger by a factor of 
4-/T 2 compared to the definition of Ref. 11911 . Note also that 
the expression for intermediate T regime, i.e. Eq (TT2b . is 
valid only in the BCS limit, and that C has a maximum at 
a particular T, since Cq, a-y, a 2 and 0:3 are all positive con- 
stants with respect to T. Using Eqs. (ITTb-dT3l> in Eq. ( fTUb . 
we find dE/da = h 2 (C - 3aiT 4 )/(47rMa 2 ), dE/da = 
h 2 (C - a 2 T 2 )/{A^Ma 2 ) and dE/da = h 2 C/{2nMa 2 ) for 
the low, intermediate and high T regimes, respectively. 

Note also that taking the derivative of the Virial theorem 
given in Eq. d29l with respect to a, and using the 1ST for the 
energy given in Eq. ( [Tol l, we obtain 



dE, 



da WirMa 2 V da dT 



(14) 



which corresponds to the 1ST for the trapping energy. Here, 



1 / 1 N 



(15) 



is the effective trapping energy, which reduces to the trapping 
energy (U) in the case of harmonic trapping potentials. 



B. 1ST for the entropy and specific heat 

The IHFT can also be used to find other thermodynamic re- 
lations. For instance, we obtain an entropy relation by taking 
the derivative of 5* = (E + F)/T with respect to A, and using 
Eqs. d2J an d ©, leading to 



dS__ d fdH\ 
~~df \~dX/ 



(16) 



In addition, we obtain a specific heat relation (at constant vol- 
ume) from C v = dE/dT = TdS/dT, leading to 



dCy = d 2 /8H\ 
d\ dT 2 \dX / 



(17) 



In Eqs. (fT&b and ( TTTb the derivatives with respect to A are 
evaluated at fixed T. We hope that these relations could 
be tested with thermodynamic measurements in strongly- 
interacting Fermi gases 12711 . 

For ultra-cold quantum gases described by the Hamiltonian 
given in Eq. (0, we can use Eq. (0 in Eqs. ( fTol l and ( TP7I ). 
leading to 
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dCy 

da 



dC 



AnMa 2 dT ' 
h 2 T 8 2 C 
i^Ma 2 df 2 ' 



(18) 
(19) 



Here, the derivatives with respect to a are evaluated at fixed 
T, and therefore, Eqs. ( fT8l and ( fT9] > correspond to the 1ST for 
the entropy and specific heat, respectively. Equation ( fT8l was 
first derived in Ref. Ql. Using Eqs. (fTTTl-dTSl in Eqs. COD 
and (O, we find dS/da = -aiT 3 / (nMa 2 ) and dC v /da = 
-3aiT 3 /(7rMa 2 ) for the low T, dS/da = -a 2 T j{2vMa 2 ) 
and dCy /da = -a 2 T/(2irMa 2 ) for the intermediate T, and 
dS/da = C/(4:TTMa 2 T) and dC v /da = -C/{2TrMa 2 T) 
for the high T regimes. Therefore, since ot\ and a 2 are posi- 
tive constants, while dS/da is negative for low T, it becomes 
positive at high T, indicating that dS/ da vanishes at a partic- 
ular temperature above Tp. Note also that dS/ da oc dCy / da 
in these regimes, and that the T dependence of C could be de- 
termined by the variation of either S or CV with respect to 



C. 1ST for the pressure 

For the Hamiltonian given in Eq. (f5]l, but in the absence of 
the potential term, i.e. a homogenous system where U(r) = 0, 
the IHFT can be used to derive the 1ST for the pressure. 

In general, via a dimensional analysis, the Free 
energy F can be written as F(r]±, . . . ,r] r ) = 
(h 2 X 2 /M)f(Xr]i, . . . , Xrj r ) where A has the dimension 
of the inverse of a length, / is a dimensionless function 
of its parameters, and r/j labels r parameters (all with the 
dimension of a length) that F may depend on for a given 
H. For homogenous ultra-cold quantum gases, since F is a 
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function of T, V, N a and a, dimensional analysis [6] requires 
that F must satisfy 



q=l 



da 



dT 



dV 



Using F = E -TS,S = -dF/dT, P = -dF/dV, and the 
1ST for the free energy given in Eq. ©, we obtain 



p _2E h 2 C 



(21) 



Therefore, the universal pressure relation in a canonical en- 
semble is of the same form as the zero temperature one HHH]. 

Taking the derivative of Eq. (l21~b with respect to a, and us- 
ing the 1ST for the energy given in Eq. ( [Tol l, we obtain 



9P_ ^ ( C[a 9C dC 
da 12nMa 2 V \ da dT 



(22) 



which corresponds to the 1ST for the pressure. We re- 
call that Eq. (l22l is derived for a homogenous system, and 
therefore, it is incorrect to compare it with Eq. ( fT4l . and 
conclude that 3VdP/da — AdE tr /da for a trapped sys- 
tem. Using Eqs. (fTTTl -lfDt in Eq. $2%, we find dP/da = 
H 2 {3C - 5aiT 4 )/(127rMa 2 U) in the BCS and dP/da = 
-bh 2 a l T A /^,TrMa 2 V) in the BEC limit for the low T 
regime, dP/da = h 2 C Q /(ATTMa 2 V) in the BCS limit for the 
intermediate T regime, and dP/da = ti 2 C / '(47t M a 2 V) for 
the high T regime. Therefore, since a.\ is a positive constant, 
while dPj da is positive in the BCS limit, it becomes negative 
in the BEC limit in the low T regime, indicating that dPj da 
vanishes at a particular scattering length around unitarity. 



D. 1ST for the atomic compressibility 

The isothermal atomic compressibility is defined as kt — 
—V~ 1 dV/dP, and it can be obtained by taking the derivative 
of Eq. (f2TT > with respect to V at constant T. Using Eq. (O 
to relate dC/dV to P, P = -dF/dV, and E = F + TS 
together with the Maxwell relation dS/dV\ T = dP/dT\ v , 
we obtain 



3 dP dP 

— =5P + a — -2T— . 
kt da dl 



(23) 



This equation also follows from a dimensional analysis of the 
pressure P which is a function of T, V, N a and a, similar to 
the analysis that lead to Eq. d20l ). The dPj da term is given by 
Eq. (1221 . and the dPj dT term can be obtained by taking the 
derivative of Eq. (l2Tb with respect to T, leading to 



dC 



W dP 



2 dT 8nMa 2 dT ' 



(24) 



where Cy — dE/dT is the specific heat at constant V. Note 
that this equation relates the specific heat to the pressure and 



contact parameter. Using Eqs. d22l and ( fT9l ) in Eq. d23l . we 
obtain 



5P 
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C 
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AT 



dC 
dT 



(25) 

which relates the isothermal atomic compressibility to the 
pressure, specific heat and contact parameter. When C = 0, 
Eq. d25l ) is satisfied for the ideal Fermi gases. The 1ST for the 
atomic compressibility can be easily obtained by taking the 
derivative of Eq. d25l l with respect to a, and using Eqs. ( l22b 
and (O. 

Note that since the compressibility is related to the den- 
sity fluctuations via the fluctuation-dissipation theorem, kt = 
FT _1 ((iV 2 ) - (N) 2 )/(N) 2 , where N is the density operator, 
it can be used to extract some thermodynamic information in 
atomic systems. Although this was proposed as early as in 
2005 B23I1 . it has recently been possible to extract this informa- 
tion for two-component Fermi gases, by measuring the den- 
sity fluctuations and atomic compressibility 124, l25ll . Com- 
bining Eq. ( fZSb with the fluctuation-dissipation theorem, pro- 
vides yet another universal relation that can be verified with 
atomic systems. Note also that the isoentropic (or adiabatic) 
compressibility k$ is related to the kt via ks/ k t = CV /Cp , 
where Cp is the specific heat at constant P. Since the spe- 
cific heats are also related to each other via Cp — Cy + 
TV n T {dP / dT) 2 , using Eq. (El, we obtain 



1 
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Kt 



TV 
~Cv~ 



dC 



2Cv 



W UwMaV dT 



(26) 



which relates k$ to Kt, Cy and C. In atomic systems, it 
is easier to measure k$ than kt, and Eq. d2"o*i i can be used to 
extract the temperarure dependence of C, given that kt can be 
extracted from the fluctuation-dissipation theorem B24l 12511 . 

Having derived the 1ST for the free, average and trap- 
ping energies, and for the entropy, specific heat, pressure and 
atomic compressibility in a canonical ensemble, next we de- 
rive the Virial theorem. 



IV. VIRIAL THEOREM FOR TRAPPED SYSTEMS 

In this section, we use the IHFT given in Eq. ([3]) to derive 
the Virial theorem J3, H2, |2^, |27l] in a canonical ensemble. 
This can be most easily achieved following the recent work 
on the zero-temperature case j26Tl . 



A. Virial theorem in a canonical ensemble 

For this purpose, consider a general Hamiltonian H = 
K + I + U, that describes N particles with arbitrary statis- 
tics in arbitrary dimensions, where K is the kinetic energy, 
/ is the interaction, and U is an arbitrary external potential. 
For ultra-cold quantum gases, the external potential is simply 
U = X^ili Ui(r{), where f7j(ri) has approximately harmonic 
dependence on the position n of the particles. 
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In general, via a dimensional analysis, U can be writ- 
ten as U(ti, rj>j) = (h 2 \ 2 /M)u(Xti, Arjv), where 
A has the dimension of the inverse of a length, and u is 
a dimensionless function of its parameters. Therefore, we 
can use the IHFT given in Eq. (01 to obtain XdE/dX — 
4(1 — Td/dT)E tr , where E tr is the effective trapping en- 
ergy defined in Eq. (fT~5T >. In addition, we can write the en- 
ergy E, via again a dimensional analysis, as E(£x, . . . , £ p ) — 
(H 2 X 2 /M)e(X£i, . . . , Xlp), where £ q labels p parameters (all 
with the dimension of a length) that E may depend on for a 
given H, and e is a dimensionless function of its parameters. 
Evaluating the derivative with respect to A for fixed values of 
£ q , we obtain XdE/dX = 2E + £ q dE/d£ q . Then, the 
Virial theorem is obtained by combining these two analysis, 
leading to 



E = 2\\-T 



d_ 

df 



E,. 



(27) 



Compared to the zero temperature expression 112611 . the main 
differences here are an extra Td/dT term in front of the po- 
tential terms, and an extra £ q term associated with the temper- 
ature. 



B. Trapped quantum gases 

In particular, for the ultra-cold quantum gases, which are 
well described by the s-wave scattering length a, Eq. d27l i re- 
duces to (1 - Td/dT)(E - 2E tr ) = -{a/2)dE/da. Note 
that for finite k c , i.e. nonzero interaction range, there would 
be an additional -{k c /2)dE / dk c = -h 2 Ck c /{4:ir 2 M) term 
on the right hand side of this equation. Furthermore, using 
the 1ST for the average energy given in Eq. (TTOb for the last 
term, solution of the resultant differential equation 12811 can 
be written as 



E = 2E t . 



h 2 c 

SnMa 



(28) 



where n is a real constant independent of T. This is the most 
general form of the Virial theorem in a canonical ensemble. 
Compared to the zero temperature expression ||4|, |5l |26J], the 
main difference here is an extra kT term. In the unitarity a — > 
±oo limit, it was shown via a dimensional analysis that the 
Virial theorem in a canonical ensemble does not have the last 
term B27I1 . and hence, we know that k vanishes in this limit, 
i.e. K a ^±oc = 0. We suspect n = for all a, however, since 
a nonzero k is allowed in general, this possibility needs to be 
clarified via other means. 

Similar to the dimensional analysis that lead to Eq. (|20b , it 
can be shown that Eq. d28"i i follows from a dimensional analy- 
sis of the average energy E supplied with Eq. ( TTOb . Note that 
E is a function of T, N a , a and the trapping frequency uj for 
a trapped system. However, applying a dimensional analysis 



to the free energy F of trapped systems, which is also a func- 
tion of T, N a , a and ui, and using Eq. (0, we obtain Eq. d28b 
with k = [29]. Therefore, we conclude that k = for all 
parameter space, i.e. 



E = 2Ef r — 



h 2 c 

S,nMa ' 



(29) 



and that the Virial theorem in a canonical ensemble is of the 
same form as the zero temperature one 10, [U, Hfl . We recall 
that taking the derivative of Eq. d29l ) with respect to a, and 
using the 1ST for the energy given in Eq. ( fTOb , we obtained 
the 1ST for the trapping energy given in Eq. ( fT4l . 



V. CONCLUSIONS 

To conclude, first we derived the isothermal Hellmann- 
Feynman theorem that is suitable for mixed states in a canon- 
ical ensemble. Then, we obtained the isothermal magnetic- 
field sweep theorems for the free, average and trapping 
energies, and for the entropy, specific heat, pressure and 
atomic compressibility of strongly-correlated ultra-cold quan- 
tum gases. We applied the sweep theorems to two-component 
Fermi gases in the weakly interacting BCS and BEC limits, 
and showed that the temperature dependence of the contact 
parameter could be determined by the variation of either the 
entropy or specific heat with respect to the scattering length. 
We also obtained the Virial theorem in a canonical ensemble, 
and discussed its implications for quantum gases. 

One of the major challenges for the experiments with ultra- 
cold quantum gases is the lack of a precise thermometry, 
and even the measure of the temperature itself for strongly- 
interacting Fermi gases is a challenging problem 13011 . On 
one hand, this makes it more difficult to perform an isother- 
mal (constant temperature) magnetic-field sweep at ultra-cold 
temperatures, compared to an adiabatic (constant entropy) 
sweep that is routinely performed in atomic systems, while 
tuning the scattering length. On the other hand, it is theo- 
retically more easier to calculate thermodynamic quantitites 
at constant temperature, e.g. the calculation of isothermal 
atomic compressibility is much easier than the isoentropic 
atomic compressibility in the BCS-BEC crossover. Therefore, 
the isothermal Hellmann-Feynman and sweep theorems that 
we discussed in this paper are probably most useful for other 
theoretical or numerical studies, and possibly for some special 
experiments where the thermometry is not an issue. 
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